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ABSTRACT 

An understanding of fuel atomization and va- 
porization behavior at superheat conditions is iden- 
tified to be a topic of importance in the design of 
modern supersonic engines. As a part of the NASA’s 
fundamental aeronautics program (the supersonics 
(SUP) project office initiative on high altitude emis- 
sions), we have undertaken an effort to assess the 
accuracy of various existing CFD models used in the 
modeling of superheated sprays. As a part of this in- 
vestigation, we have completed the implementation of 
a solution procedure for the modeling of superheated 
sprays into the national combustion code (NCC), 
and then applied it to investigate the following three 
cases: (1) the validation of a flashing jet generated 
by the sudden release of pressurized R134A from a 
cylindrical nozzle, (2) the differences between two su- 
perheat vaporization models were studied based on 
the predicted results of both hot and cold flow cal- 
culations involving a Parker-Hannifin pressure swirl 
atomizer, (3) the spray characteristics generated by 
a single-element LDI (Lean Direct Injector) combus- 
tion experiment were studied to investigate the dif- 
ferences between superheat and non-superheat con- 
ditions. For the validation case, the predicted re- 
sults provided reasonable agreement for both gas and 
droplet velocities but underpredicted the droplet sizes 
by a considerable measure in the central region. For 
the second case, the superheat vaporization models 
of both [12-14] and [15] produced similar results for a 
hot flow. For the third case, both the superheat and 
non-superheat calculations produced mostly similar 
results for gas temperature and velocities but the su- 


perheat calculation produced substantially faster va- 
porization as expected. 

NOMENCLATURE 

Bk Spalding mass transfer number 

B t Spalding heat transfer number 

C p specific heat, J/(kg K) 

d droplet diameter, m 

dt time increment, s 

h specific enthalpy, J/kg 

k thermal conductivity, J/(ms K) 

Ik mixture latent heat of evaporation, J/kg 

lk,eff effective latent heat of evaporation, 

J/kg (defined in Eq. (14)) 
m liquid mass flow rate, kg/s 

rhk, flash droplet vaporization rate under 

flash evaporating conditions, kg/s 
rhk,t droplet vaporization rate due 

to heat transfer, kg/s 

n.fc number of droplets in kth group 

Nil Nusselt number 

P pressure, N/m 2 

Pr Prandtl number 

P 3 at saturation pressure, N/m 2 

rfe droplet radius, m 

R u gas constant, J/(kg K) 

Re Reynolds number 

Sh Sherwood number 

Sfc droplet radius-squared ( = r|), m 2 

Sc Schmidt number 

SMD Sauter mean diameter, m 

t time, s 

T temperature, K 
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Tf, boiling temperature, K 

7i. kth droplet temperature, K 

U gas or liquid velocity, m/s 
x axial distance, m 

y radial distance, m 

A p pressure drop in the injector, N/m 2 
p dynamic viscosity, kg/ms 

v kinematic viscosity, m 2 /s 

p density, kg/m 3 

9 spray cone angle, deg. 

Subscripts 

/ fuel 

g gas 

inj injector 

l liquid 

r radial coordinate 

s droplet surface 

t time 

x axial or x-coordinate 

y y-coordinate 

z z-coordinate 

Superscripts 

mean, or average 
flow rate 

INTRODUCTION 

Our previous work on spray modeling can be 
found in Refs. [1-10]. The spray solution procedure 
facilitates multi-dimensional spray/gaseous combus- 
tion calculations on massively parallel computers and 
unstructured grids. It also allows calculations to be 
performed in conjunction with the scalar Monte Carlo 
PDF (Probability Density Function) method [2] . The 
objective of our present study is to extend our solu- 
tion procedure to the modeling of superheated sprays 
because of its importance in the design of modern su- 
personic engines, scramjet, and ramjet afterburners 
[11-15]. It is because higher heat-load requirements 
may lead to the use of the the same fuel as a coolant 
and also sometimes the nozzles are required to oper- 
ate at low back pressures. Also, under normal gas- 
turbine operating conditions, it is estimated that a 
small fraction of the liquid fuel may be flash boiled 
[12]. There are also some reported cases of flash re- 
lated engine performance problems in gasoline direct- 
injection internal combustion engines [12]. Although 
flash evaporation is considered to be detrimental to 


engine performance under normal circumstances, it 
can have some potential benefits: it is known to 

produce a fine spray with enhanced atomization, in- 
crease effective spray cone angle, and decrease spray 
penetration [11]. Because of its importance in some 
NASA-related applications, we have undertaken an 
assessment study to establish baseline accuracy of 
various existing CFD models used in the modeling of 
a superheated spray. As a part of this effort, we have 
implemented a solution procedure into NCC based on 
the modeling approach developed by [12-15]. More 
specifically, we have implemented a superheat vapor- 
ization model and other modifications pertinent to 
the atomization modeling developed by [12-14]. Also, 
we have completed the integration of a computer 
module that was developed by UTRC (United Tech- 
nologies Research Center) and UMASS (University of 
Massachusetts) [15]. The UTRC /UMASS computer 
module contained a separate superheat vaporization 
model, a flash-induced atomization breakup-model, 
and a calculation approach that was developed to de- 
termine the JP8 two-phase properties [15]. In or- 
der to establish a baseline accuracy of the solution 
procedure, we have undertaken a validation involving 
a flashing jet generated by a cylindrical nozzle [16]. 
Also, we have performed some calculations involving 
a Parker-Hannifin pressure swirl atomizer in order 
to evaluate the differences between the superheat va- 
porization models of [12-14] and [15], and a single 
element LDI reacting case in order to evaluate the 
differences between the superheat and non-superheat 
conditions [16]. 

In this paper, we first describe complete details 
of the solution procedure implemented into NCC fol- 
lowed by a presentation of the results. 

ATOMIZATION MODELING 

The atomization refers to a process of the liq- 
uid jet breakup into droplets. There are many pro- 
cesses associated with the liquid jet breakup. They 
can broadly be classified into two regimes: (1) the in- 
ner nozzle flow, and (2) the external nozzle flow. In 
the inner nozzle flow, several factors (such as injector 
type, geometry, and size) influence the conditions at 
the injector exit (such as the velocity, the initial sheet 
or jet thickness, and the angle of droplet dispersion). 
Right now, we rely on widely-used correlations in de- 
termining the initial spray conditions at the injector 
exit. For a better description, a more accurate anal- 
ysis is needed which takes into account the physics 
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associated with inside bubble growth, cavity forma- 
tion, and internal turbulence. In the external nozzle 
flow, the liquid becomes unstable under the influence 
of aerodynamic instabilities and finally breakup into 
droplets. When the maximum amplitude of the most 
unstable wave reaches a critical value, some liquid is 
stripped of the main liquid core in the form of pri- 
mary droplets. These droplets may further breakup 
into smaller droplets due to a process known as the 
secondary droplet breakup mechanism. 

Among the primary atomization models, the 
following are widely used: (1) the sheet breakup pri- 
mary atomization model which was developed by [17- 
18] has application in the modeling of sprays gen- 
erated by pressure swirl atomizers, (2) the blob jet 
primary atomization model which was developed by 
Ritz et al [19-21] has application with plane orifice 
nozzles, and (3) the BLS (Boundary-Layer Stripping) 
primary atomization model which was developed by 
Khosla et al [22] has application for a liquid jet in 
a cross flow. Among the secondary droplet breakup 
models, the following are widely used: (1) the TAB 
(Taylor Analogy Breakup) secondary droplet breakup 
model which was developed by O’Rouke et al [23] 
has application for lenticular-shaped droplet defor- 
mations, and (2) the ETAB (Enhanced TAB) sec- 
ondary droplet breakup model which was devel- 
oped by Tanner [24] has application for high-pressure 
diesel-engine flows. 

Most of the coding required for both the pri- 
mary atomization and secondary droplet breakup 
models was provided by CFDRC. These models are 
applicable under both non-superheat and superheat 
conditions. In the superheated regime, the thermo- 
physical liquid properties are replaced by the two- 
phase properties of a superheated fluid [15]. 

Flash Atomization Modeling 

There are three components to the atomization 
modeling of a superheat spray: (1) primary atom- 
ization, (2) secondary droplet breakup, and (3) flash 
induced atomization. Further details on the primary 
atomization and secondary droplet breakup models 
can be found in [1]. For the superheated sprays, we 
have to deal with a third component known as flash 
induced atomization. We have followed the modeling 
approach of Lee et al [15] for the flash induced atom- 
ization. In what follows we describe some details of 
their modeling approach. 

Flash atomization is a result of rapid evapora- 
tion associated with strong nucleation and bubble 


growth occurring inside of a non-equilibrium two- 
phase fluid [15]. When the void fraction within the 
rapidly evaporating fluid reaches a certain critical 
void fraction, the inner liquid core would shatter into 
pieces. It seemed to have a value of 0.61 based on 
the research identified in [15]. Incidentally, this con- 
dition for criticality also coincides with the packing 
limit of spheres beyond which it can not accommo- 
date any more bubbles and results in shattering [27] . 
In order to calculate this state of critical limit, they 
employed the homogeneous relaxation model (HRM) 
for calculating the mass fraction of gas, x, also known 
as quality. 


§ = o = 6» 0 eV; (i) 

where the values for a (= —0.54), /3 (= —1.76), and 
d 0 (= 6.51.10 o4 s) are taken from [28], and the void 
fraction, e, and the degree of superheat, (j>, are defined 
as follows: 


Pi ~ P . , = Pb(h) - P 

- Pv P C — Pb{h) 


Once the critical limit for the void fraction is 
reached, the resulting droplets produced as a result 
of flash atomization are obtained by a droplet number 
density correlation developed by Kawano et al [29] . 


n = 4e (T v^ub) e at + b (3) 

From Eq. (3), the mean droplet size, R, is cal- 
culated by equating total liquid mass before and after 
the droplets are formed. 


R 3 


1 - e e 
(4/3)7rn 


(4) 


The solution for Eqs. (1) to (4) as well as the 
JP8 fluid properties needed in this calculation were 
obtained by the computer coding provided by UTR.C. 
Also, it is noteworthy that the properties for JP8 were 
predicted by means of a surrogate model of Violi. 
Further details can be found in [15]. 

The flash induced atomization is combined with 
the primary and secondary droplet breakup models 
following the approach of Lee et al [15]. The aero- 
dynamic instability mode will atomize the liquid core 
when the surface disturbance reaches a certain critical 
value and the flash induced atomization is responsi- 
ble for shattering the liquid core when the void frac- 
tion reaches the critical value. Therefore, the control 
mechanism depends on whichever mechanism is first 
responsible for reaching the critical condition. 
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Some Additional Aspects of Superheat on Primary 
Atomization 


primary atomization models needs further investiga- 
tion. 


Flash evaporation is known to produce a fine 
spray with enhanced atomization, increase effective 
spray cone angle, and decrease spray penetration [11]. 
Here, we consider the effects of flash evaporation on 
the sheet breakup primary atomization model by fol- 
lowing the approach of Zuo et al [12]. Its effect on 
the initial droplet size generated immediately after 
the first ligament breakup, d ls , is given as a function 
of both engine pressure and a superheat parameter 
as follows: 


VAPORIZATION MODELING OF A 
SUPERHEATED DROPLET 

In what follows we first provide some details 
of the two superheat vaporization models that were 
implemented into NCC. 

Superheat Vaporization Model of Zuo, Gomes, 

& Rutland [12], and Schmehl & Steelant [13-14] 


= d,„(-^)°- 27 [i-x(^) 0 ' 135 ] o < *(%M 0 ' 135 < i 

■* atm * * 

( 5 ) 

where di n is the corresponding droplet droplet size 
under normal evaporating conditions without flash 
evaporation, and y is defined as a superheat param- 
eter as follows: 


I{T k ) - I(T b ) 

TO 


( 6 ) 


where / is the internal energy of the liquid. Its value 
varies between 0 < x < 1 with y = 0 referring to zero 
flash evaporation and \ = 1 to full flash evaporation. 

In Eq. (5) , the increase in d. is due to an increase 
in engine pressure by a factor of (py— ) 0 ' 27 is based 
on an experimental correlation obtained from Lefeb- 
vre [30] . It reflects the influence of chamber pressure 
on wave propagation as it damps wave growth. But 
the decrease by a factor of (1 — ) 0135 ) j s due 

to a significant reduction in droplet size caused by 
both cavitation and bubble growth under flash evap- 
oration conditions. It was introduced based on the 
experimental data obtained from VanDerWege et al 
[31] and Reitz [32]. 

As the liquid approaches boiling, it also causes a 
substantial decrease in both intact liquid core length 
and core droplet size leading to a modification of the 
nominal cone angle, 9 , as given by 


9 = 9 n + (144 — 9 n )x 2 (7) 

where 9 n is in degrees for a spray vaporizing under 
normal conditions without flash evaporation. This 
correlation was developed by based on the experimen- 
tal data of Reitz [32], Both these correlations were 
developed in conjunction with the sheet breakup pri- 
mary atomization model but their validity with other 


It is based on an extension of the classical D 2 - 
theory. In the classical evaporation model, the ther- 
mal energy needed for evaporation is mostly furnished 
by the external heat transfer from the surrounding 
gas. Under superheat conditions, the characteristic 
vaporization time resulting from the external heat 
transfer from the surrounding gas is of the same order 
of magnitude as that resulting from the flash evap- 
oration. The energy needed for vaporization at the 
droplet surface is partly provided by the superheat 
energy stored within the droplet and it is controlled 
by the droplet internal heat transfer. The modeling 
approach differs from the classical droplet vaporiza- 
tion models in three important ways: (1) under su- 
perheat conditions, the droplet surface mass fraction, 
Yf s , approaches unity as the droplet surface temper- 
ature is maintained at the corresponding liquid boil- 
ing temperature; (2) under superheat conditions, all 
the external heat transfer from the surrounding gas 
is made available to the vaporization process with no 
apparent increase in the droplet surface temperature; 
and (3) the flow of fuel vapor imparted by flash va- 
porization partly counterbalances the flow generated 
by external heat transfer and may significantly reduce 
the energy transferred from the surrounding gas. 

Based on the governing equations of conserva- 
tion for an isolated spherically symmetric droplet, 
Zuo et al [12] and Schmehl and Steelant [13-14] 
showed that the total evaporation rate, m k , can be 
calculated as 


dl k — dlk^flash 4” rhk,t (8) 

where the flash boiled vaporization rate, ihk, flash . , is 
given by 


TTlk, flash 


47T r\a s 


(T k - T b ) 
Ik 


(9) 
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where T k is the internal droplet temperature and the 
overall heat transfer coefficient, a s (= kJ/s m 2 °K) 
is given by the Adachi correlation [26] : 

= 0.76(7), - T,,) 0 ' 26 (0 < T ks -T b < 5) 

a s = 0.027(T fc - T 6 ) 2 ' 33 (5 < T ks - T b < 25) (10) 

= 13.8(7), - T 6 )°' 39 (T ks -T b > 25) 


The Adachi correlation is valid over a wide range of 
superheat conditions. 

The vaporization rate due to external heat 
transfer, m k ,t, in Eq. (8) is given by 


= 27rr fc - 


-'p 1 + 


Nu 


_ M l +( 1 + w )Bt 
: m k t 


( 11 ) 


where the Spalding heat transfer number, B t , and the 
Nusselt number, Nu, are given by 


B t 


C p (T g - T ks ) 

lk,eff 


(12) 


Nu = 2(1 + 0.3 Re 1/2 Prl /3 ) (13) 

where the effective latent heat of vaporization, h,eff , 
is given by 


h,eff = h + 47t 


Vfc 

TOfc 



(14) 


which is an useful parameter as it represents the total 
energy loss associated with the latent heat of vapor- 
ization in addition to the the heat loss to the droplet 
interior. Finally, the corresponding droplet regression 
rate, is given by 


dsk _ rhk 

dt 27t r k pi 


(15) 


This model is valid over an entire range of su- 
perheat conditions as long as there is some amount of 
superheat energy available within the droplet ( T k > 

T b ). 


Superheat Vaporization model of Lee et al [15] 


This modeling approach is also based on an ex- 
tension of the classical _D 2 -theory but it differs from 
the previous model of [12-14] in several important 
ways: (1) neglects the effects of internal nucleation 
and convection, (2) instead of invoking the Adachi 
correlation for the flash induced vaporization rate, 
the heat loss from the droplet interior is modeled by 


di = interior T S (P 00 ) /g\ ^ formulation is valid 

^ ffc,i=o ’ v ' 

only when T, x > Tf ru b ( Poo ) , which means for posi- 
tive Spalding heat transfer numbers. 

With the assumptions employed, Lee et al [15] 
arrive at the following equation for the evaporation 
rate, rh v . 

T \ = ™>vL-qi (16) 

where 

T = C P gT/L r , m v = (4irr 2 p g u)/(4Tr\ g rk,t=o/C P g), 
L = L k /L r . The solution for the above equation is 
obtained by employing the standard Newton’s itera- 
tion. 

It is noteworthy that because of the assump- 
tions invoked in this model, the applicability of this 
model may be limited in several important ways: 
(1) our experience shows that the assumption of 
Too A Tbub(Poo) becomes too restrictive in most 
non-reactive calculations and also in some regions of 
most reactive spray calculations, and (2) the qi term 
is modeled based on the assumption that the overall 
droplet life is very short (sub millisecond range) but 
in most moderate superheat conditions the droplet 
lifetimes are on par with those vaporizing under non- 
superheat conditions. So its use is likely to be re- 
stricted for the modeling of smaller droplets produced 
after flash induced atomization. 

The solution for Eq. (16) and all the JP8 fluid 
properties needed in this calculation were also ob- 
tained by the computer coding provided by UTRC. 

Vaporization Model Valid Under Non-Superheat 
Conditions 


Under moderate initial superheat conditions, 
only a fraction of the vaporization takes place under 
superheat conditions ( T k > Tb) and the remainder 
takes place under non-superheat evaporating condi- 
tions (Tfc < Tb). So there is a need to revert back to 
a vaporization model valid under stable evaporating 
conditions when the internal droplet temperature ap- 
proaches the boiling temperature. In the present cal- 
culations, the vaporization rate under normal evapo- 
rating conditions is evaluated by means of a simplified 
classical Z? 2 -theory: 


m fe = 27 t r k p g D fgs Sh ln( 1 + B k ) (17) 

where the Spalding mass transfer number, B k , is 
given by 
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Bk 


( Vfs ~ Vf) 
(! - Vfs) 


(18) 


where yy s is the mass fraction of the fuel species at 
the the droplet interface, yf is the mass fraction of the 
fuel species in the surrounding gas, and the Sherwood 
number, Sh, is given by 


Sh = 2(1 + 0.3 Re l/2 Sc 1 J 3 ) (19) 

Internal Droplet Temperature Calculation 


Our experience with the validation studies 
showed us that there is a definite need to include 
a calculation involving the internal droplet tempera- 
ture valid under both superheat and normal evapo- 
rating conditions. In our present calculations, it was 
evaluated by means of a simple infinite conductivity 
model. 


dT k _ 3 [lk,eff ~ h\ dsk 
dt 2C p i r l dt 

if T k < T b , and 


(20) 


dT k 

dt 


3cr s 

TkPlCpi 


C T k - T b ) 


( for the superheat model of [12 — 14]) 


or 


dT k _ 3A; dT 

dt r k piC p i dr 


( for the superheat model of [15]) 


if T k > T b 


the sudden release of pressurized R.134A from cylin- 
drical nozzles of different sizes. In our present calcu- 
lations, we have focused mainly on a case involving 
the 1mm nozzle for which the most extensive data 
were reported. 

It is based on a 2d-axisymmetric geometry: 2d- 
axisymmetric. All the calculations were performed on 
a grid with 79101 triangular elements (non-reacting, 
turbulent air flow). 

At the start of jet dispersion, the pressurized 
liquid would undergo some expansion due to the en- 
trainment of surrounding air and vaporization of the 
pressurized liquid depending on the degree of super- 
heat. Because of this expansion, the effective size 
of the expanded jet would be larger than the orifice 
exit. Based on the experimental observation, the ini- 
tial size of the jet at the point of jet breakup was 
estimated to be about 9 mm. The initial gas velocity 
was prescribed based on the experimental correlation 
reported by Yildiz et al [34] . The initial gas temper- 
ature 298 K, turbulence kinetic energy, k = 1 & = 
10 m 2 /s 2 . The liquid is R134A, a refrigerant. The 
injector has a diameter of 0.001 m and is located in 
the middle of a center-body configuration. The initial 
droplet velocity was estimated to be 34 m/s based on 
the pressure drop across the injector. The initial liq- 
uid temperature is 293 deg. K which corresponds to 
a superheat temperature of 46 deg. C. 

The results for this case were reported in our 
earlier paper [16]. The predicted profiles for both gas 
and droplet velocities show a reasonable agreement 
with the measured data, and exhibit a self-similar 
behavior as reported in the literature. The predicted 
droplet sizes mainly fall short of experimental data 
in the central region of the spray. Away from the 
central region, the predicted sizes are in reasonable 
agreement with the experimental behavior. Both the 
predictions and experimental data showed that a con- 
siderable amount of vaporization takes place under 
non-superheat conditions. 


RESULTS & DISCUSSION 


Parker-Hannifin Pressure Swirl Atomization 


A Flashing Jet Generated by a Cylindrical Nozzle 

This investigation is under taken mainly to test 
the following models: (1) the superheat vaporization 
model of [12-14], (2) the blob-jet primary atomization 
model [19-21], & (3) the ETAB secondary droplet 
breakup model [24-25]. The experimental data for 
this validation was provided by Yildiz et al [33-34]. 
It contained data for two-phase jet flows generated by 


This investigation is undertaken mainly to test 
the following models: (1) Sheet-breakup primary 

atomization model [17-18], & (2) ETAB secondary 
droplet breakup model [24-25]. The experimental 
data for this validation was provided by Ref. [35]. 
It is based on a 2d-axisymmetric geometry. For this 
case both the grid and the initial conditions were pro- 
vided by CFDRC (non-reacting, turbulent air flow). 
The experimental conditions for the gas phase are 
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specified as follows: (1) the axial velocity is assumed 
to be uniform at 1.2 m/s, (2) the radial velocity pro- 
file is prescribed by means of an interpolation from 
the following three points: y = 0.019 m, v = -0.9m/s, 
y = 0.1 m, v = -0.5 m/s, & y = 0.3 m, v=0, (3) the 
gas temperature = 300 K, & (4) the turbulence ki- 
netic energy, k = 1 & = 10 (m 2 /s 2 ). The liquid phase 
conditions are specified as follows: (1) the liquid fuel 
is C12H23, (2) the injector has a diameter of 0.0012 m 
and is located in the middle of a center-body config- 
uration, (3) the liquid mass flow rate is 0.00125 kg/s, 
(4) the internal pressure drop across the atomizer is 
344500 Pa, & (5) the initial liquid temperature is 300 
K. 

The experimental conditions were based on the 
non-superheat evaporating conditions. The reported 
experimental results involved droplet axial velocities, 
droplet radial velocities, and droplet sizes. The com- 
parisons for this case were also reported in our earlier 
paper [16]. The predicted results were in reasonable 
agreement but mostly underestimate the the experi- 
mental data. 

We have also completed two more calculations 
for the same case of Parker Hannifin but performed 
under non-reacting superheat conditions. In one, the 
SVM (Superheat Vaporization Model) of [12-14] was 
employed and in the second the SVM of [15] was 
employed. In both calculations, the properties of 
C12H23 were replaced by the JP8 properties of a su- 
perheated fluid. They were calculated based on the 
coding provided by UTRC [15]. The specified initial 
liquid temperature was raised to 520 deg. K. For the 
JP8, the boiling temperature at the normal pressure 
is about 468 deg. K. 

The global structure of the non-reacting spray 
obtained from both the two calculations are shown in 
Fig. 1. The differences between the two calculations 
are considerable: (1) The results from the SVM of 
[12-14] showed a temperature variation ranging from 
270 to about 430 deg. K. Initially there is an increase 
in the gas temperature near the spray region from 
the superheated droplet. While some of the smaller 
droplets vaporize completely under superheat condi- 
tions, some of the larger droplets cool off below the 
boiling temperature after depleting all of the super- 
heat energy from the droplet interior. Further va- 
porization from those larger droplets tend to lower 
the temperatures in the surrounding gas below the 
corresponding droplet surface temperature. This ef- 
fect can be seen downstream in a small central region 
where temperatures fall below the initial specified gas 


temperature of 300 deg. K. This behavior of evapo- 
rative cooling is similar to what we have seen in the 
validation case of Yildiz et al [16]. (2) The results 
from the SVM of [15] show the gas temperatures to 
range between 300 to 355 deg. K. This calculation 
also shows a slight increase in the gas temperature 
in the central region near nozzle exit. In this cal- 
culation the gas temperatures mostly stay above the 
initial specified gas temperature of 300 deg. K. It is 
because this superheat model is only valid for posi- 
tive Spalding numbers; therefore, makes its applica- 
bility more difficult for the modeling of a non-reacting 
spray. 

We have also completed two more calculations 
for the same case but performed under reacting super- 
heat conditions. Again in one calculation the SVM 
of [12-14] was employed and in the second the SVM 
of [15] was employed. All other conditions are main- 
tained the same as in the non-reacting superheat cal- 
culations. 

The global features of the reacting spray ob- 
tained from the reacting calculations are shown in 
Fig. 2. The spray particle distribution appears to 
be similar in both calculations. The gas temperature 
ranges between 300 to 2100 deg. K. In both calcula- 
tions most of the combustion takes place downstream 
of the spray in a predominantly diffusion-controlled 
flame. However, in the calculation involving Lee et 
al [15], the high temperature region extends further 
upstream into the inner central region of the hollow- 
cone spray. 

Swirl Stabilized Reacting Spray Generated 
by a Single Element LDI Experiment 


This investigation is undertaken to test the fol- 
lowing models: (1) the Sheet breakup primary at- 
omization model [17-18], & (2) the TAB secondary 
droplet breakup model [23]. The experimental data 
for this validation was provided by Ref. [36]. Fig. 
3 shows the test section of the experiment reported 
in [36]. The air passes through a six helicoidal vane 
swirler followed by a converging-diverging venturi be- 
fore entering the combustor. The air mass flow rate 
is 0.49 kg/min., and at the inlet, air speed=20.4m/s, 
T=294.3K, P=101325Pa. The fuel nozzle is located 
at the middle of the centerbody and it is a Parker- 
Hannifin, 90-deg. hollow cone, pressure swirler. The 
liquid fuel is Jet-A, the liquid mass flow rate 0.025 
kg/min, and the initial liquid temperature 475 K. The 
internal pressure drop across the atomizer is 350000 
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Pa which leads to an initial liquid velocity of about 
22 m/s. 

The results are based a 3D calculation with a 
grid size of 833024 elements. It makes use of the 
following models for the gas-phase: the standard 

k-epsilon turbulence model, a a single-step reduced 
chemical kinetic mechanism, and the eddy-breakup 
turbulence/chemistry model. The experimental con- 
ditions were based on the non-superheat evaporat- 
ing conditions. The reported experimental results in- 
volved gas and droplet velocities (axial, radial, and 
tangential), & droplet sizes. The comparisons for this 
case were also reported in our earlier paper [16]. For 
this case, the predicted gas-phase results underesti- 
mate both size and location of the central recircula- 
tion region but the predictions provide better agree- 
ment for both droplet size and velocities. 

We have also completed two more calculations 
for this case involving superheat conditions. In one, 
the SVM of [12-14] was employed, and in the second, 
the SVM of [15] was employed. In both calculations, 
the liquid properties of the jet-a fuel were replaced by 
the two-phase properties of JP8. Also, the specified 
initial liquid temperature was raised to 520 deg. K. 

Figs. 4 to 10 show the results for both gas 
and droplet velocities, and droplet sizes. Shown here 
are the comparisons between the experimental data, 
the predicted obtained from the corresponding non- 
superheat experimental conditions, as well as the pre- 
dicted results obtained from the two different super- 
heat calculations. The comparisons for gas axial ve- 
locities at different axial locations are shown in Fig. 
4, and the comparisons for the radial and tangential 
velocities are shown in Figs. 5 and 6, respectively. 
The gas velocities obtained from the superheat con- 
ditions are similar to those obtained from the non- 
superheat conditions and they all underestimate both 
size and location of the central recirculation. The 
predictions could be improved by taking into account 
of the unsteady nature of fluid flow and also by in- 
corporating a better chemical kinetic mechanism and 
chemistry/turbulence model. 

The comparisons for the droplet sizes are shown 
in Fig. 7 at four different axial locations. The com- 
parisons between the experimental data and the pre- 
dictions from the non-superheat calculation are in 
good agreement. It is also noteworthy that the ra- 
dial spreading rate of the spray is well predicted. 
As you move progressively further downstream, more 
droplets start disappearing from the the central re- 
gion. It is because of the rise in the gas tempera- 


ture to near flame temperature which leads to faster 
evaporation. But the comparisons differ significantly 
with the results obtained from the superheat calcu- 
lations. The results from both the superheat vapor- 
ization models are mostly similar. At the first axial 
location, the predicted droplet sizes from the super- 
heat conditions are higher than the experimental data 
in the the central region. The reason for this increase 
could partly be attributed to a decrease in the liq- 
uid density caused by a rapid bubble growth in the 
superheated fluid. This would lead to a decrease in 

the liquid breakup length/time because of its . — 

y P a 

dependence which in turn leads to an increase in the 
initial product droplet size because of its inverse rela- 
tionship to the liquid breakup time. But at the next 
(second) axial location, the drop sizes mostly fall be- 
low the experimental results. Further downstream, 
the effect of superheat conditions is more evident as 
the predicted droplet sizes fall significantly below the 
experimental data. 

Fig. 8 shows the droplet mean axial velocity 
comparisons at six different axial locations. Shown 
are the predicted results as well as the experimental 
data for two different droplet-size groups: 0-15 mi- 
crons and 60-75 microns. The data for the 0-15 mi- 
crons represent the velocities of the smallest droplet 
sizes and the data for the 60-75 microns represent the 
velocities of the largest droplet sizs. Our predicted re- 
sults are expected to fall somewhere in between the 
two droplet-size groups but they are expected to fall 
closer to the 60-75 microns as the predicted results 
are based on the mean-size (SMD) droplet represen- 
tation. The comparisons with the non-superheat cal- 
culation are mostly in good agreement and fall mostly 
closer to the large-size droplet group but for some un- 
derestimation observed in the middle of the first axial 
location, and the outer region of the second axial lo- 
cation. The results from the superheat calculations 
are also mostly in similar agreement in the first three 
axial locations but the differences grow in the last ax- 
ial locations where fewer droplets are present due to 
faster vaporization under superheat condtiions. Also, 
the axial velocities are underpredicted in the last two 
axial locations. Fig. 9 shows the droplet radial ve- 
locity comparisons. There seems to be a reasonable 
agreement between the predicted and experimental 
results but for some underestimation observed in the 
middle of the first axial location, and the outer re- 
gions of the second and third axial locations. Also, 
it is interesting to note that the radial velocities for 
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the small-size droplets are much higher than the cor- 
responding larger-size droplets in the first two axial 
locations. It is because the smaller droplets tend to 
follow the gas velocity where the radial gas veloci- 
ties are much higher because of swirl-induced recir- 
culation. The mean droplet tangential velocity com- 
parisons are shown in Fig. 10. The predicted re- 
sults fall mostly closer to the experimental data of 
the large-size droplets at the first three axial loca- 
tions. They mostly underpredict the experimental 
data in the next three locations but the superheat 
results show a slightly better improvement. It is be- 
cause smaller droplets resulting from increased vapor- 
iztion due to superheat tend to follow the gas velocity. 
The comparisons could perhaps be improved by pro- 
viding an additional swirl component to the initial 
liquid velocity as no swirl component was specified in 
the present calculations. 

CONCLUDING REMARKS 

The paper describes complete details of a solu- 
tion procedure that was implemented into NCC for 
the modeling of a superheated spray. The model- 
ing approach was used in the validation of a flashing 
cylindrical jet, and in the model-accuracy assessment 
studies based on a pressure swirl atomizer to examine 
the differences between the two superheat vaporiza- 
tion models, and a single-element LDI reacting case 
to evaluate the differences between the non-superheat 
and superheat conditions. 
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(b) SVM of Lee et al [15]. 


Fig. 1 Global structure of a non- reacting spray generated by a Parker- ffannifin 
pressure-swirl atomizer under superheat conditions. 
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(b) SVM of Lee et al [15]. 

Fig. 2 Global structure of a spray flame generated by a Parker-Hannifin 
pressure-swirl atomizer under superheat conditions. 
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Fig. 3 Test setup of the single-element LDI experiment 
(taken from Cai et al [36]). 
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Fig. 4 Gas mean axial velocity distributions. 
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Fig. 5 Gas mean radial velocity distributions. 
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Fig. 6 Gas mean tangential velocity distributions. 
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Fig. 8 Spray mean axial velocity distributions. 
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Fig. 9 Spray mean radial velocity distributions. 
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Fig. 10 Spray mean tangential velocity distributions. 
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